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Abstract 

The classical approximation provides a non-perturbative approach 
to time-dependent problems in finite temperature field theory. We 
study the divergences in hot classical field theory perturbatively. At 
one-loop, we show that the linear divergences are completely deter- 
mined by the classical equivalent of the hard thermal loops in hot 
quantum field theories, and that logarithmic divergences are absent. 
To deal with higher-loop diagrams, we present a general argument 
that the superficial degree of divergence of classical vertex functions 
decreases by one with each additional loop: one-loop contributions 
are superficially linearly divergent, two-loop contributions are super- 
ficially logarithmically divergent, and three- and higher- loop contri- 
butions are superficially finite. We verify this for two-loop SU(A^) 
self-energy diagrams in Feynman and Coulomb gauges. We argue 
that hot, classical scalar field theory may be completely renormalized 
by local (mass) counterterms, and discuss renormalization of SU(A^) 
gauge theories. 
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1 Introduction 



The classical approximation 111] is a useful tool for the study of infrared properties 
of quantum fields at high temperature 0, ^ ^, |7|, ^, |^ , which may be applied 
to calculate non-perturbative phenomena such as the Chern-Simons diffusion rate 
[p!0| , \n\ , |T^] (relevant for theories of baryogenesis JlSl) and the dynamics of the 



electroweak phase transition [0], as well as real-time (plasmon) properties of hot 



non-abelian gauge theories The classical theory is expected to be a good 

approximation at low-energy because the classical limit ^ ^ and the low-energy 
limit of the Bose- Einstein distribution function n yield the same result: 

1 1 

nM = T-rr — 7 — = nciM, huJk<t:T, (1.1) 



where (-^k = vk^ is the frequency at wave-number k and (3 = 1/T the inverse 
temperature. Classical correlation functions are determined by a set of field 
equations in Minkowski space and a thermal average over the initial fields at some 
(arbitrary) initial time. In perturbation theory, classical vertex functions can be 
obtained by taking the limit — *■ of the quantum expressions, which amounts to 
the replacement (|1 . 1|) of the Bose-Einstein distribution function by the classical 
distribution function. The resulting ^'s in the denominator are compensated by 
a positive power of h's arising from the loop counting of the diagrams under 
consideration, such that in the classical limit a non-trivial expression (and loops) 
remain. 

The replacement (|1 . 1|) is a good approximation for infrared-dominated di- 
agrams, but it changes the ultraviolet behavior of the theory and introduces 
classical (Rayleigh- Jeans-type) divergences. When the classical theory is con- 
sidered as a low-energy effective theory, these divergences can be regularized by 
introducing a cut-off of the order of the temperature, A ~ T/^. Since in a weakly 
coupled theory the temperature is large compared to dynamically generated en- 
ergy scales such as g^T, the resulting cut-off dependences are a direct reflection 
of the divergences of the classical theory. The general strategy to improve the 
effective theory is to include counterterms that reduce the cut-off dependence. In 
particular, if a complete set of counterterms can be specified, the cut-off may be 
sent to infinity and the theory is renormalized. It is clear that a knowledge of 
the divergences is necessary to determine the appropriate counterterms. We will 
assume that these divergences will be tractable in perturbation theory. 

In the case of a scalar field theory the divergences have been studied in 
classical perturbation theory for the two-point function up to two loops and the 
four-point function up to one- loop [§], 0, |]. It was found that the one-loop resp. 
two-loop correction to the self-energy is linearly resp. logarithmically divergent, 
and that the one-loop correction to the four-point function is finite [|^. In 3 + 1 
dimensional gauge theories on the other hand, the attention has mainly been 
restricted to the classical equivalent of the quantum hard thermal loop (HTL) 
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expressions |]T6|, 0, |T8| , which introduce hnear divergences in the classical theory 
P, ^. Numerical studies using a HTL improved effective theory P, |T9|, ^ can 
be found in 22, 23 1. An analysis of the divergences in the classical theory 



that goes beyond the HTL limit at one-loop, or to higher loops, has not yet 
been performed for gauge theories. Our aim in this paper is therefore to give a 
more complete analysis of the divergence structure of hot, real-time classical field 
theory. 

A different kind of effectively classical theory beyond the HTL regime was 
constructed in ||2^, by integrating out the scales T and gT in a leading log ap- 



proximation. It takes the form of a Langevin equation and is free from ultraviolet 
divergences [^]. For a numerical implementation, see p6|.Q However, our focus 
in this paper is on classical Yang-Mills theory as it stands, without any further 
integrating out to construct an effective theory. 

We shall argue that both in SU(A^) gauge theory and in scalar field theory with 
0^ and 0^ interaction terms the divergences are restricted to one- and two-loop 
(sub) diagrams. It will be shown that classical one-loop diagrams that correspond 
to HTL's in the quantum theory lead to linear divergences, while other one- loop 
diagrams are finite in the classical theory. Also we present a general argument 
that two-loop diagrams can at most give logarithmic divergences. This is explic- 
itly verified for two- loop self-energy corrections in SU(iV) and scalar theories. 

The paper is organized as follows. Classical one-loop diagrams are analyzed in 
the next section, and diagrams with two loops and more in section |^. In section 
1^ the possibility of absorbing the divergences with counterterms is discussed. 
The final section contains the conclusions. Throughout the paper expressions 
for classical diagrams are obtained by taking the h ^ limit in the quantum 
expressions. The validity of this is discussed in more detail in appendix 0, where 
a set of classical Feynman rules is presented for hot scalar fields. 



2 One-loop 

2.1 Linear divergences: classical HTL's 

The one-loop linear divergences of the classical theory are closely related to the 
(quantum) hard thermal loops discovered by Braaten and Pisarski [|1^] (see also 
[|17| , ^). For instance, the divergent part of the classical self-energy in SU(A^) 
gauge theory can be obtained as the classical limit of the HTL self-energy |^, ^. 

^This effective theory has been rederived with the help of classical transport theory, using 
the concept of classical colored point particles ||2^ . It should be clear that we study classical 
fields instead of classical particles. 
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To be specific, the spatial part of tlie retarded HTL self-energy readsQ 

^i^TL.AP) = -26^VnN [ ^ kkn'M ^ ^\ , (2.1) 

where here and in the following the external frequency p° is taken real with a 
small imaginary part to obtain the retarded self-energy, i.e. p*^ = Re(p°) +ie, and 

As usual in the HTL approximation, the radial and angular integration decouple 
and the radial integration determines the plasmon frequency 

^pi = -^/^^X dken'{k) = -g'N—. (2.3) 

The classical self-energy corresponding to (|2.1| ) is obtained by taking the h — * 
limit, before the integration over k is performed. This simply amounts to 
replacing the Bose-Einstein distribution function by the classical distribution 
function, as in ( |1 . 1| ) . The classical self-energy is non- vanishing, since the h in 
the prefactor of ( p.l|) is compensated by the h in the denominator of the classical 
distribution function. The resulting radial integral is linearly divergent and to 
handle this we introduce a cut-off in the classical distribution function nci(co'k) — > 
nc\{uJk)6{A — k). This particular way of introducing a momentum cut-off in loop 
integrals does not lead to problems with gauge invariance, which can be most 
easily understood from the gauge propagator of Landshoff and Rebhan |^0| and 
is explained in appendix [A.2| . The result is a linearly divergent classical plasmon 
frequency 
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2 



^pi,ci - ir^g'NTA. (2.4) 

The relation between the quantum plasmon frequency ( |2.3| ) and the classical ana- 
logue (|2.4|) is that the Bose-Einstein distribution function effectively introduces 
a cut-off of the order of the temperature on the integration, A ~ T/h. Since 
the angular integration is completely decoupled, the dependence on the external 
momenta of the linearly divergent contribution to the classical self-energy and 
HTL self-energy are equal.0 All of this is well-known [0, 

Hard thermal loops are the leading contributions to vertex functions for soft 
external momenta ~ gT. Power counting reveals that one-loop diagrams, 

with any number of external gauge fields, contain a HTL contribution. The fact 
that the external momenta are small compared to the internal momentum ~ T 

^Loop momenta will generically be denoted with K = (fc°,k) and external momenta with 
P = (P°,P)- Furthermore, K'^ = -/cq + k'^, k = |k| = Wk, and k = k/k. 

•^At least with a (perturbative) continuum-like regularization as employed here. On a spatial 
lattice, this is not the case [p|, p|, |3l[ . 
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allows for several simplifications in the calculation of HTL's. As a result all HTL's 
are proportional to the plasmon frequency squared ( p.3|) |TF 



Divergences in classical field theories have a similar behavior, since here also 
the internal momenta ~ A are much larger than the external momenta. In fact, 
all classical HTL's have the proportionality factor (|2.4| ). Therefore, all classical 
HTL's are linearly divergent. 

Other one-loop contributions in the quantum theory are smaller by a factor 
p/fc ~ p/T. In the classical limit these subleading contributions give a factor 
p/k ^ p/ which reduces the degree of divergence. Therefore we may conclude 
that all linear divergences at one-loop are given by the classical HTL's. 

2.2 No logarithmic divergences 

Next we will argue that there are no logarithmic divergences at one-loop in the 
classical theory. Firstly, we discuss one particular example in SU(A") gauge theory 
explicitly, which is the spatial part of the self-energy in the Feynman gauge. A 
convenient starting point is the expression in the quantum theory, which reads 



Hg'(P) = 5''^g^hN j 



d^k \ 2n{uji,) + 1 Aij 



1 1 



-^[^(cuk) - n(Lc;p+k)] 



+ t^k + Wp+k - CUk - t^p+k 

1 1 



with 



- C^k + t^p+k P^ + UJi,- CJp+k 

(2.5) 

1 



Skikj + hpikj + ?,kiPj + A{p - pQ)gij - 2piPj . (2.6) 



This diagram contains of course the HTL self-energy ( ^.1| ). As before, the classical 
expression is obtained by taking h to zero. The non-thermal contribution from 
the "1" in the first and second line vanishes as h goes to zero. 

From the previous section we know that contributions to the self-energy ( |2.5| ) 
are at most linearly divergent. The classical limit of the momentum- independent 
tadpole-like contribution in the first line is indeed linearly divergent. For the 
contribution proportional to Aij, it implies that the contributions bilinear in the 
external momenta, i.e. the terms proportional to piPj or p'^Sij, can only give ul- 
traviolet finite contributions, and that the terms linear in the external momenta 
(terms proportional to kipj or Pikj) may give logarithmic divergences. The con- 
tributions proportional to kikj may contain logarithmic divergences besides the 
linearly divergent contributions as well. 

To obtain the linearly and logarithmically divergent contributions we expand 
the integrand in 1/k, so that we can estimate the ultraviolet behavior of the 
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integrand by power counting. The contribution from the second hne reads 



A. 



4ci;kt^p+k 



1 



^ ^ f n„(A:) + (p ■ k) jn'M ' i^n^iik) + 0{k 



Ak"^ [k 



k 



k^ 



(2.7) 



The first term on the second hne, with Aij oc kikj, is part of the HTL contribution. 
The second term between curly brackets, and the first term with Aij oc Pikj, kiPj, 
contain the contributions proportional to k~^, and these may give a logarithmic 
divergence after integration. However, it turns out that these contributions are 
odd under the transformation k —>■ — k and therefore they vanish upon integra- 
tion. The other terms, including those indicated with 0{k~^), are ultraviolet 
finite by power counting. 

Similarly, the third line can be expanded, and after some algebra it can be 
written as 



Ai 



4a;ku;p+k 



n(uj 



p+k 



)] 



1 



1 




- cuk + Wp+k + u;k - t^p+k 



p-k)<i(A:) + -(p-k)ViW 



-^<,{k) 



2P^ 



(P • k)^ 



- (P • k)^ 



+ 0{k 



-4\ 



(2i 



The first term on the second line, again with Aij oc k^kj, is part of the HTL 
contribution, and is proportional to The other terms contain a contribution 
proportional to k~^, which after integration could yield a logarithmic divergence. 
However, just as in the previous case these contributions are odd under the trans- 
formation k —>■ — k and they vanish upon integration. The remaining terms are 
ultraviolet finite. 

Therefore, we conclude that there is no logarithmic divergence in the spatial 
part of the retarded classical self-energy in the Feynman gauge. In a similar 
manner, we have also verified that the spatial part of the three-point vertex 
contains no logarithmic divergences. 

The reason for the vanishing of possible logarithmically divergent contribu- 
tions lies in the behavior of the self-energy and the vertex functions under parity 
(P) and time reversal (T). The spatial part of the self-energy discussed here is 
invariant under p — p, and p^ —>■ —p^ in combination with complex conjuga- 



tion (i.e. p + ie —{p + i^) in (|2.5| )). The point is that the expansion in 1/k 



turns out to be an expansion in PT odd (dimensionless) functions of p° and p. 
Since the linearly divergent HTL contributions to the self-energy are even under 
P and T, the logarithmically divergent contributions are odd and should therefore 
vanish. This argument extends to the temporal part of the self-energy as well as 
to other vertex functions. 
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Finally we would like to remark that the vanishing of logarithmic divergences 
holds in general Coulomb or covariant gauges, since the corresponding gauge 
fixing term does not break PT invariance, and the same argument can be applied. 



2.3 Classical self-energy: explicit result 

The analysis presented above is useful for a general understanding. However, 
in some cases it is possible to actually calculate the loop integrals and avoid an 
expansion in 1/k. Here we give one of those explicit results in SU(A^) theory. 

We calculate the diagonal {ii) part of the classical one-loop retarded self- 
energy in the Feynman gauge in appendix 0, and the result reads 



nti(^) 



1^-?^ ^ IT 



3p2 



-I m — 



P 



— p 47r \ " 2p p'' — p 

The real and imaginary parts can be obtained in the usual way, using 



(2.9) 



In 



pO _|_ p 
pO _ p 



In 



p^ + p 



pO _ p 



-mO{p^ -pI) 



(2.10) 



The linear divergence is precisely the equivalent of the hard thermal loop 
contribution, which follows from the replacement TA/vr^ T'^/{6h). The finite 
terms are exactly equal to the terms linear in T that are obtained in a high 
temperature expansion in the quantum theory, as can be checked explicitly 
3^.0 There are no other terms. The p° — limit equals the well-known result 
from the quantum theory in the Feynman gauge |^ 



nti(o,p) 



3pT 



(2.11) 



Note that in this limit the leading-order (gauge-dependent) behavior is completely 
determined by classical physics. 

To conclude the one-loop analysis, the above described situation can be un- 
derstood also directly by keeping h in the high-temperature expression of the 



quantum theory. The high-temperature expansion then has the form |£9 



^H_i(P) + THo(P) + (h In ^ j niog(P) + hU, (P) + ' 



(2.12) 

where fi is the renormalization scale. The term proportional to is the HTL 
part, which turns into the linearly divergent term when h ^ 0, and the second 
term in this expansion is the finite term in the classical theory. All the other 
terms vanish when h ^ 0. 



'Up to some typographic errors. 
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3 Two-loop and beyond 



3.1 Degree of divergence 

In this section we study the degree of divergence of higher-loop diagrams in the 
classical theory. In the first part we shall argue that the superficial degree of 
divergence of the self-energy decreases by one with each loop, starting with the 
one-loop linear divergence. Then we will check this statement explicitly for a 
number of diagrams. We shall argue that the same is true for classical vertex 



functions in section 3.3 



To make the argument for the self-energy, we start with the following basic 
assumption: in the high-temperature limit the retarded self-energy in the quan- 
tum theory scales according to its dimension, i.e., the quantum retarded gluon 
self-energy behaves as 

n^,(p) = r%,(p7p,p,<7) + r2o(p/r), (3.1) 

for high temperatures, fixed external momentum and frequency, and a renormal- 
ization scale of the order of the temperature /i ~ T. This assumption consists 
of two parts: The contribution of diagrams with hard momenta 7^ ~ T on all 
internal lines gives a contribution to the self-energy. Contributions that are 
excluded in (|3.1|) are of the form g'^^T'^{T/ P)"^ for m > and with L indicating 
the number of loops. For fixed external momenta and high temperatures such 
terms become larger than the one-loop (HTL) contribution g'^T'^, so they invali- 
date a loop expansion. Therefore the assumed absence of these contributions can 
be re-expressed by saying that we assume that hard modes are perturbative. The 
other part of the assumption is that also diagrams with soft internal momenta 
give a contribution. This relies on the belief that infrared divergences are 
controlled by induced masses, which are proportional to the temperature, such 
as the electric and magnetic masses in SU(A^) gauge theories. 

Let us then consider a classical contribution to the self-energy containing M 
distribution functions. To be able to compare the degree of divergence of such 
a contribution with the quantum expression, we regard the temperature in the 
quantum self-energy as a particular ultraviolet cut-off, and using the assumption 



( |3.1|) we count the degree of divergence as 2. Since every classical distribution 
function gives rise to an extra energy in the denominator when compared to the 
quantum diagram]^ the classical contribution to the self-energy with M distribu- 
tion functions has then a superficial degree of divergence 2 — M. 

To complete the argument, we now use that the number of distribution func- 



tions M can be related to the number of loops L in the following manner p4| , p5 



One way to obtain the retarded self-energy is by using the imaginary-time or 

^In the ultraviolet regime of a loop integral the quantum (Bose) distribution function can 
be approximated as exp —fihuj and acts as a cut-off function. On the other hand, the classical 
distribution function remains proportional to l/uj. 
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Matsubara formalism |T8[ . One first performs the sums over the discrete loop 
frequencies and then analytically continues the external frequencies to real values 
with a small positive part to incorporate the appropriate retarded boundary con- 
ditions. In the imaginary-time formalism the number of loops equals the number 
of Matsubara frequency summations. Using the method of contour integration 
to perform these sums, each sum gives rise to one 'coth' function, either with 
positive or negative energy. Explicitly, each sum gives a factor ^| 



1 shoo 11, 

2 coth^ = n(s^) + - = + -], s = ±. (3.2) 

Hence, the resulting expressions are of the form of spatial momentum integrals 
over Bose-Einstein distribution functions, where the number of distribution func- 
tions is equal to or less than the numbers of loops. The classical limit can now 
be taken by replacing n{uj) + ^ — > T/{hu), such that the h's counting the loops 
cancel against the from the distribution functions. After taking the classical 
limit, only the leading term, which has as many distribution functions as loops, 
remains and the number of classical distribution functions M in a given diagram 
is counted by the number of loops, M = L. Note that this applies not only to 
the self-energy diagrams but to vertex functions as well. It follows then that the 
superficial degree of divergence of a classical diagram is given by 2 — L, such that 
the classical one-loop contribution to the self-energy is superficially linearly di- 
vergent, the two-loop contribution is superficially logarithmically divergent, and 
higher-loop contributions are superficially finite. 

3.2 Two-loop self-energy diagrams 

We now want to verify the general argument of the previous section for the two- 
loop self-energy diagrams appearing in SU(A^) and scalar field theory. We do not 
discuss diagrams which have a one-loop self-energy subdiagram (and hence also 
a linear sub divergence), but we concentrate on the two-loop diagrams as shown 
in fig. m. Furthermore, since we are only interested in the structure of ultraviolet 
divergences, i.e. in power counting, we do not need to make a distinction between 
gauge field propagators in the Feynman gauge and ghost propagators in the loops. 

Let US, clS db first relatively simple example, take the two-loop setting-sun 
contribution (a) to the retarded self-energy as it appears in A^'^-theory (with 
A = g"^) and SU(A^) gauge theory. It reads 

u^^Hp) = -iQ'h? I ^ y - 

^ ^ 6^^ W (27r)3 7 (27r)3 23cUkCUk'C^ki P° + scuk + + sic^ki 
{ [1 + ^(scjk)] [1 + n{s'uj]^')] [1 + n(siu;ki)] - n(scJk)?^(sW)^(sit^ki)}, (3.3) 
where tUki = cUp-k-k', and the sum is over all s's being ±. 
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(a) (b) (c) 

Figure 1: Two-loop diagrams. The setting sun diagram (a) and diagram (b) are 
discussed in section pT2| , and diagram (c) is treated in appendix |C[ 

Note that the product of three distribution functions drops out. It is then 
clear that the classical limit of fB.3|), 



ii,, [n - J J X. 



X 



(27r)3 J (27r)3 23^k^k'^ki 

1 ( , r2 \ ^ 

Si + s' + s , (3.4) 



contains products of two classical distribution functions, in accordance with the 
statement that the number of loops equals the number of distribution functions. 
We now estimate the degree of divergence by power counting and take the loop 
momenta k, k' ~ A. The integral measures give two contributions ~ A^, and 
all single energy denominators 1/uj give a factor 1/A. The energy denominator 
that contains p° will produce, for generic large loop momenta k, k', a hard energy 
denominator ~ 1/A. It can only produce a soft energy denominator when there 
is a cancellation, which is in the special case that k ~ ±k', depending on the 
signs of s, s' and si [l^ . However, for these special configurations the integral 
over phase space is restricted so that this will not alter the degree of divergence. 



We will use this estimate for energy denominators with three hard energies |]16 
below as well. 

By power counting we therefore establish that this contribution is logarith- 
mically divergent, as expected. This is also the result obtained in [|, |^, |[, where 
the classical setting sun diagram was analyzed in detail and it was shown that 
in fact the logarithmic divergence can be separated and is independent of the 
external momentum and frequency. 

It should be noted that the setting sun diagram (as well as the diagrams dis- 
cussed below) contains an infrared divergence for vanishing external momentum 



. For massless A(/>^ theory, this can be cured by resumming the effective ther- 
mal mass, arising from the one-loop tadpole diagram. This has only an effect on 
the soft infrared modes, and does not interfere with the ultraviolet behavior of 
the classical diagram we investigated above. 

The next example we treat is the two-loop diagram (b) in fig. ||, which appears 
in SU(A^) and in scalar 0'^-theory. This particular diagram is more delicate. 
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and it is instructive to carry out the procedure described above in detail. We 
will verify explicitly that in SU(iV) theory (the spatial part of) this diagram is 
logarithmically divergent in the Feynman gauge. 

Since we are only interested in the degree of divergence of the diagram, we 
may ignore the color and Lorentz structure of the diagram. To indicate the 
momentum-dependence of the four vertices in the gauge theory, we will insert a 
factor (A;)^^-. The precise form of the momentum insertions is unimportant for the 
power counting. 

We have found it convenient to calculate this diagram in the imaginary-time 
formalism, and after performing the sums over the Matsubara frequencies, the 
diagram can be written as 



! (- 

-\- SUJ -\- S2OJ2 I - 



SS'SIS2S3 

X 



pO + S'CU' + SsUJs p'^ + sou + S2U2 I —S^U^ -\- S2U2 + SlUi 

{\ji{siUJi) l][n(s2a;2) l]n(s3a;3) - n(sia;i)n(s2a;2)[n(s3a;3) 1 
1 



X 



p° S3a;3 + suj — siUJi 
([n(s3a;3) l\[n{suj) l\n{siUJi) - n{s2,oj^)n{suj)[n{siUJi) 1]) 
1 



X 



P° -|- S'uj' -\- S2OJ2 — SiUJi 

{\ji{s'uj') l][n(s2<^2) + '^\n{siUJi) — n{s' uj')n{s20J2)[n{siUJi) 1]) 
1 



X 



S'u' — SU -\- S\U\ 

{\n{siUJi) -\- l\[n{s'uj') -\- l\n{suj) — n{siUJi)n{s' uj')\n{suj) -\- 1])!? 

(3.5) 

where we have used the shorthand notation 

and the sum is over all sign factors s = ±1. Again the products of three dis- 
tribution functions drop out. The corresponding classical integral n|^j,j(P) may 
be obtained by taking the h ^ Q limit, which amounts to neglecting the con- 
stants and single distribution functions and replacing all distribution functions 
that appear in products of two by classical distribution functions. 

We will now consider the large k,k' ^ A behavior of the classical diagram as 
we did for the setting sun diagram, by looking at the various factors in n^^|,[(P) 
and naively combine those to obtain an indication for its degree of divergence. 
First of all, each integration measure contributes (Pk ~ and the factor {k)fj 
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is proportional to A^. Each of the energies in the denominator on the first hne 
gives a contribution Xjuj ~ 1/A, such that this factor leads to a contribution 
1/A^. Each classical distribution function gives a factor 1/A as well. 

The other energy denominators require a bit more care. All energy denomina- 
tors between the curly brackets contain three large energies that will generically 
not cancel, as in the case of the setting sun diagram. These therefore contribute 
with a factor 1/A. The two energy denominators in the second line may produce 
a 'soft' energy denominator for specific combinations of the sign factors, namely 
for S3 = — s', and S2 = —s. For example, the first denominator may give 

A°, (3.7) 



+ s'(c^k' - t^p-k') + s'k ■ p 

similar to what happens in the one-loop case. This gives us four possibili- 
ties: both energy denominators produce a soft contribution, only one of them 
is soft and the other is hard, or both are hard. Putting all these estimates to- 
gether, we obtain in the first case, with two soft denominators, the naive result 
nSi,ci(^) ~ A6A4A-5A-2A-1 ~ A2, which is a quadratic divergence. With one 
soft denominator we find n^^),j(P) ~ A, a linear divergence, and with two hard 

contributions n|^|.j(P) ~ A°, the expected logarithmic behavior. However, from 
the general argument we expect a logarithmic divergence only. 

The reason for this mismatch is that this naive power counting doesn't treat 
the distribution functions correctly. In the one-loop (HTL) case, often differences 
of statistical factors appear. In the classical theory, these lead to a different 
ultraviolet behavior and hence change the power counting. Therefore we take a 
closer look at the two-loop diagram to see whether a similar thing occurs here 
as well. We denote the (naively) quadratically divergent piece, with S3 = — s' 
and S2 = — s, with n-^|,j(P). To re-estimate the divergence, we put the external 
momentum in the energy denominator with three large loop-energies (i.e. in the 
second, fourth and sixth line of ( |3.5| )) equal to zero, since for generic large k, k' 
the denominator does not vanish.]^ Taking the external momentum equal to zero 
can in fact be seen as the zeroth order term in an expansion in the external 
momentum. The first order term, linear in the external momentum, is treated 
in appendix y. The naively quadratically divergent contribution can now be 
written, after flipping s\ to —S\ in the term on the sixth line, as 

(b) (T3\ ^1 2f-\2 /" f /m4 '^1 -'■ 



(27r)3 J (27r)3 ^ 2^luu;'uj,u;2UJ3 + s'{uj' - ^3) 



P° + s{uj — LU2) S'oj' — SUJ + SiOJ\ 

(3. 



^Again, the region where it does vanish is only a restricted part of phase space and is 
excluded in the argument for power counting. 
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We redo the power counting for Ilij'^i{P). The thing to notice is that indeed two 
differences of two classical distribution functions have appeared, and for hard 
loop-momenta 

[nci(st^p_k) - ridisujk)] ~ -s(k ■ p)?2^i(u;k) ~ A"^ (3.9) 

Both differences give one extra power of 1/A, compared to the naive power 
counting employed before. The conclusion is therefore that Ilf!l^^i{P), instead 
of being quadratically divergent, is only superficially logarithmically divergent, 
as expected by the general argument. 

Note that the classical limit of diagram (b) may contain a linear divergence 
from a HTL (three-point) subdiagram. The linear divergence occurs, e.g. in 
contribution ( p.8| ), whenever {k)fj ~ k^k' or {k)fj ~ kk'^. However, at this stage 
we are not interested in divergences caused by one-loop subdiagrams since we 
study only the superficial degree of divergence. 

Potentially, there are also superficial linear divergences in the classical limit 
of (p.5|). These are worked out in appendix 0. In this appendix we also discuss 
the other self-energy contribution (c), which is naively linearly divergent as well. 
It turns out that they are all in fact logarithmically divergent, in accordance with 
the general argument of the preceding section. 



3.3 Higher-order vertex functions 

We now extend the argument to general vertex functions. At zero-temperature 
we know that the degree of divergence of a Feynman diagram decreases with the 
number of external lines. In a real-time classical theory at non-zero temperature 
this is not the case. We already saw that the linear divergences at one-loop occur 
for diagrams with any number of external gauge field lines. Therefore we do not 
expect that the two-loop contributions to three- or higher-point functions are 
finite. 

To argue what happens for vertex functions with more loops, we use the 
real-time Feynman rules which are presented for scalar field theory in appendix 
A.l|. We employ Feynman rules in which two type of propagators appear, the 



temperature-independent retarded propagator Gq and the thermal two-point 
function 5*0 that contains the thermal distribution. It is useful to recall here 
their explicit representation 

G^{K) = E ^T^r--^- , So{K) = Y: n^^{su^)^6{k'-suj^). (3.10) 

Starting from the classical retarded self-energy with L loops (and hence M = L 
thermal propagators), generalized retarded ra-point functions with L loops can be 
obtained by adding retarded Green functions in the loops, using the vertices (a) 
and (c) shown in fig. B of appendix |A.1|. Note that thermal propagators cannot 
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(a) (b) 




Figure 2: Two-loop diagrams in the real-time formulation that contribute in the 
classical limit. Full lines are thermal propagators and dashed-full lines retarded 
propagators. 



be added in the loops, since then the number of distribution functions M is no 
longer equal to the number of loops, which is required by the argument given in 
section and is needed to have the cancellation of h in the classical diagrams. 
Note that this also implies that all integrals over the zeroth components of the 
loop momenta can trivially be performed with the help of the on-shell delta 
functions in the thermal propagators. To continue, in the case of a gauge theory, 
every additional (momentum-dependent) three-point vertex gives an additional 
factor K (we do not need to be more specific for the power counting argument 
presented below). Hence the total effect of adding one external line using a 
three-point vertex is an additional factor K times a retarded propagator 

. (3.11) 

From the viewpoint of power counting, the first factor is of order A°, and the 
second factor can be of order A'^ or 1/A, depending on whether a soft or hard 
energy denominator results, after the integrals over the on-shell delta functions 
in the thermal propagators have been performed. 

This leads us to give the following general argument: in the case that the 
propagator in ( p.llj ) is soft, the additional external hne will not change the degree 
of divergence, compared to the diagram without the additional line. On the other 
hand, when the energy denominator turns out to be hard, when the extra vertex 
is a 4-point vertex, or in scalar field theory, where the momentum K in the 
numerator is absent, additional lines will always lower the degree of divergence. 
Using the result for the two-point function, this implies that higher-point vertex 
functions are superficially logarithmically divergent by power counting (at two- 
loop) or finite (at higher- loop). 
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Figure 3: Two-loop diagram with three external lines. 



There is one slight complication in this general argument. In the self-energy 
considered in the previous section, the logarithmic divergence was the result of a 
subtle cancellation between quadratically (and linearly) divergent contributions. 
The question is whether this subtle cancellation is not spoiled by adding an 
external line. Although a complete analysis of two-loop vertex functions is beyond 
the scope of this paper, we will check explicitly in one particular case that the 
cancellation indeed still occurs. 

This analysis can be done most conveniently using the real-time Feynman rules 
of appendix ^ We start by presenting in fig. ^the classical two-loop contribution 
to the self-energy (b) in the real-time formalism. The integral over the zeroth 
components of the loop momenta can easily be performed using the on-shell delta 
functions in the thermal propagators, and we have verified that this yields indeed 
the classical limit of ( p.5| ), which was calculated in the imaginary-time formalism, 
as expected. 

We want to add one external line to obtain a diagram as in fig. ^. In the case of 
the self-energy that we discussed in the previous section we found that the naively 
quadratically divergent contribution (|3.8|) does not contain a distribution function 
at energy uji = cuk-k'- That means that in terms of the real-time diagrams no 
diagram with a thermal propagator on the line shared by the two loops contributes 
to ( p.8| ). Hence we do not need to consider the addition of extra lines to the third 
and fourth diagram. Let's now see how an additional three-point vertex of type 
(a) in fig. |^ and an additional retarded Green function can be added to the first 
two diagrams in fig. ^. It turns out that for each diagram (a) and (b) there are 14 
possibilities to do this. A closer look, however, reveals that not all diagrams are 
needed to establish a cancellation of the naive quadratic divergence. For example, 
a combination of the two diagrams that are shown in fig. ^ is sufficient to obtain 
a difference between distribution functions that reduces the degree of divergence 
to a logarithmic one. 

Indeed, the sum of the most divergent part of the diagrams in fig. ^ yields 



^ J (27r)3 J (27r)3 ^"'^^'^ 2^uiulu^^^, 



SS'Sl 

1 1 1 



pI + spi-k + sp2 ■ k gO - s'q ■ k' SicUk-k' - - 
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Figure 4: Two-loop contributions to the classical 3-point vertex function in the 
real-time formalism that combined yield a logarithmic degree of divergence. 



[raci(su;pi-k) - riciisujk)] [nci(s'u;q_k') - nci(s'a;k')] , (3.12) 

with Q = P1 + P2. The factor (A;)fjfc has been included to account for the momen- 
tum insertions from the vertices in a SU(A^) gauge theory, and the factor h'^ arises 
from loop-counting. We had to expand also the single energy denominators, such 
as l/cijp^-k, in external momenta. Compared to the self-energy expression ( |3.8| ) 
the vertex function has one extra factor (|3.11|) with a soft energy denominator as 
anticipated. After power counting, taking into account ( |3.9| ), we may conclude 
that in this particular combination the addition of one external line does not spoil 
the reduction from a quadratic divergence to a logarithmic divergence. 

It will be interesting to make explicit checks for other three- (and higher) 
point vertex functions with two loops as well, but without a clever method to 
combine the different contributions this seems to be out of the question. 



3.4 Other gauges 



To verify the general argument in section ^Tl| that two-loop diagrams are log- 
arithmically divergent, we have estimated in sections |3.2| and p.3| the degree of 
divergence of some two-loop diagrams in the Feynman gauge. Here we want to 
argue that the estimates in the Feynman gauge extend to general Coulomb gauges 



le 



The gauge propagator in a general Coulomb gauge with gauge parameter ac 
reads 

1 1 K Ki, 

^t^u = -^T^M + '^MO'^z.Op + ^4 (3.13) 

with the transverse projector Tjj(k) = 6ij — kikj /k"^ ,Tqq = Tqi = Tjo = 0. 

First we realize that the external momentum dependence in the transverse 
projector may be neglected T^iy(p — k) ~ T'/ii/(k) when the integration momentum 
k is large. In the power counting of a diagram we may estimate T^^, ~ 1, and 
we see that a diagram with all transverse propagators has the same degree of 
divergence as the same diagram in the Feynman gauge. Since the 00-component 
and the gauge dependent part of the propagator cannot give a soft denominator 
like (|3.7|), we can also neglect the external momenta in these components, they 
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are then estimated as Therefore diagrams containing these components of 
the propagator will not have a larger degree of divergence. We conclude that the 
degree of divergence of a certain diagram is the same in the Feynman gauge and 
in a general Coulomb gauge. We stress that this does not necessarily imply that 
the logarithmically divergent contribution is gauge independent as is the case for 
the linear divergences, this remains a subject for further study. 

Finally we like to remark that in general covariant gauges it is not expected 
that individual diagrams obey the power counting of section 3J-, but rather the 
sum of the diagrams with a certain number of loops. 



4 Renormalization 
4.1 Scalar field theory 

The general analysis given in the previous sections yields the following result for 
a scalar field theory with 

^nt = / d'x + iA0^) (4.1) 

as interaction term: time-dependent classical scalar field at finite temperature is 
renormalizable. 

To see this in more detail, let's start at one loop: first of all, the tadpole 
diagram, which is the only HTL in the quantum theory, is linearly divergent, 
as expected. But since this divergence is trivially independent of the external 
momentum, it can be canceled with a mass counterterm p, ^. The self-energy 
correction with two 0^ vertices is finite. This can be seen in a number of ways: 
the general analysis revealed that in gauge theories such a diagram is linearly 
divergent. In a scalar theory, the vertices are momentum independent and this 
brings down the superficial degree of divergence with 2, which makes the diagram 
is finite. In other words, it is not a HTL. It follows also from a explicit calculation 
0, where the one- loop correction to the vertex function in a A0^-theory was 
studied, which is equal to the one-loop self-energy we are discussing now. It 
was found there that the classical result is the leading order term in a high 
temperature expansion in the quantum theory. Finally, in the quantum theory 
one-loop contributions to higher n-point functions are no HTL's either, therefore 
these diagrams will be ultraviolet finite in the classical limit. The reason is 
that the extra retarded Green functions that appear in the loop bring down the 
superficial degree of divergence with at least one, and this cannot be compensated 
by momentum- dependent vertices as in the gauge theory. 

At two loop, the setting sun contribution to the self-energy is logarithmically 
divergent as expected from the general analysis. This has been calculated explic- 
itly in 0, H . However, it has been shown there that the divergent part can be 
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isolated and is independent of the external momentum and frequency: again the 
divergence can be absorbed with a mass counterterm 0. For the other contribu- 
tions to the self-energy that were discussed in detail in the previous section, we 
note the following: in the gauge theory, diagram (b) contains a power counting 
factor [kY ~ due to the momentum-dependent vertices. In the scalar case 
this is of course absent. Hence the naive quadratically divergent contribution in 
the gauge theory is superficially finite in the scalar case, even without any need 
for subtle cancellations. The same is true for diagram (c), here a power counting 
{ky ~ is absent and the classical diagram is immediately superficially finite 
as well. We conclude that for the two-loop self-energy corrections (neglecting di- 
agrams that contain self-energy subdiagrams) the setting sun diagram is the only 
diagram that is superficially divergent. The absence of momentum-dependent 
vertices also simplifies the analysis of n-point vertex functions. It was argued in 
section that the superficial degree of divergence goes down at least by one if 
more retarded Green functions are added in the loops in a scalar theory. This 
can be immediately applied here. The complication in gauge theories, i.e. the re- 
quirement of a subtle cancellation, is not needed for the vertex functions either, 
because even without the cancellations, the diagrams are superficially finite. 

Since the general analysis already shows that three loops and higher are finite, 
we arrive at the conclusion formulated at the beginning of this section: classical 
scalar field theory at finite temperature can be renormalized with merely a mass 
renormalization. For more details on this and on e.g. the choice of finite parts, 
see [§. 

We would like to stress that the analysis has been limited to simple correlation 
functions, as discussed above. Other quantities, such as transport coefficients like 



the viscosity |]38[, which require expectation values of composite operators [B9 



are not covered here. 



4.2 One-loop renormalization in SU(A/^) theory 

In section § we have seen that the divergences at one-loop are given by the classical 
HTL's. Below we will present counterterms for these one- loop divergences. 

Since we are interested in a classical theory it is appropriate to consider the 
equations of motion 

5^^S=[D,,F'"^] = 6r, (4.2) 

with the covariant derivative D^, = + igA,^, and the field strength F'^'^ = 
d'^A^ — d^A"^ + ig[A'^, A'^]. The source consist of a finite renormalization and 
a divergent part. We write 

sf = jL-jL- (4.3) 

The divergent source j^j^ has to be chosen such that it cancels the linear di- 
vergences 0. The identification of the linear divergences as classical HTL's 
then implies that the divergent source must equal the classical HTL induced 
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source: j^'j^ = Jhtl ci- We take the finite renormalization such that the clas- 
sical theory ( [4.2|) is consistent with the quantum results. This can be seen as 
the real-time analog of matching employed in static dimensional reduction and 



the construction of effective 3-dimensional theories [^]. Therefore, the finite 
renormalization should equal the quantum HTL's: jg^ = Jhtl- ^ (classical) 
perturbative expansion this renormalization prescription means that any classical 
HTL (sub) diagram is replaced by the corresponding quantum HTL diagram. 

The HTL effective action is non-local, but the introduction of an auxiliary 
field allows for a local formulation of the HTL equations of motion [0. The 
sources then read 

= fHTL = ^g'^N J -0^^n'Mv^W{x,^r), 



■fi 

Jdh 



with the velocity f ^ = (1, v) and v = k. The velocity-dependent auxiliary field 
W{x,\) satisfies the equation 

[^;^D^^y(x,v)] = v■E(x), (4.5) 

with E the (chromo-)electric field. 

From (|4.4|) we see that the two sources can be combined into one expression 

^ [n'ico^) - n',,{u;^)]v^W{x, v). (4.6) 

As usual in the HTL approximation the radial integration can be performed 
independently, which yields 

6f = 36uj^ I — t;W(x,v), (4.7) 

J 477 



with the linearly divergent coefficient |42, 20 



= - (4.8) 



(see section We see that the whole set of non-local linear divergent vertex 
functions can be renormalized by the adjustment of one parameter. This is a 
special feature of the continuum HTL's. On a lattice the counterterms are more 
complicated |P, |, since the lattice breaks rotational invariance and the ve- 



locity on the lattice is not equal to the speed of light. The effective theory with 
continuum-like counterterms (|4.7|J4.8| ) has been studied numerically in |2^. An 



important observation is that the equations ( |4.2| , |47^ , [4.8| ) have conserved 
energy [HO] , such that a thermal average over initial fields can be defined and the 
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equations are a part of a proper classical statistical theory. Furthermore, it is 
well-known that in the classical theory only the combination g'^T appears and % 
is absent. However, since cUpi oc g'^T'^/h, T/h is introduced as an independent, 
non-trivial scale in the effective theory described above. 

That no linear divergences occur in retarded vertex functions calculated with 



the classical theory ([4.2| , ^77| , |47q ), can be seen by introducing a background 
field and integrating out the classical fields. This generates in the HTL approxi- 
mation the classical induced source which is precisely what is subtracted on the 
r.h.s. of (|4.2| ). It would be nice to see in a perturbative calculation that these 
counterterms are sufficient to absorb the linear (sub) divergences in (superficial 
logarithmically divergent) two-loop diagrams, but this has not been attempted 
here. 



5 Conclusion 

Classical thermal field theories contain ultraviolet divergences. In an analysis 
of classical vertex functions, we found that at one-loop only linear divergences 
occur, which come from classical HTL's, i.e. the classical equivalences of the 
HTL's in the quantum theory. Furthermore we argued that for ra-point vertex 
functions with arbitrary n, the degree of divergence decreases with the number of 
loops. This implies that two-loop contributions are (superficially) logarithmically 
divergent and higher loops are superficially finite. This may be compared with 
static dimensional reduction, where the L-loop contribution to the self-energy has 
also a degree of divergence 2 — L. The difference is that in the static limit higher- 
point vertex functions are less divergent than the self-energy. Indeed, the static 
theory is a superrenormalizable field theory and a finite number of counterterms, 
like a one- and two-loop mass counterterm, suffices. 

The consequences of our findings are the following. Since three and higher- 
loop diagrams are superficially finite, these are infrared dominated. Therefore, 
they are in principle calculable in the classical theory. The loophole is of course 
the possible occurrence of divergences in (one or two-loop) subdiagrams. To deal 
with these divergences, counterterms have to be introduced. In the scalar case 
the divergences occur only in the self-energy and are momentum independent, 
therefore a mass renormalization is sufficient to obtain a cut-off independent the- 
ory. This may be useful for a numerical approach to time-dependent problems, 
such as the dynamics of the phase transition and/or topological defects in a (com- 
plex) scalar field theory. In SU(A^) gauge theories the divergences are momentum 
dependent, nevertheless a renormalization of the plasmon frequency ( [4.8| ) takes 
care of the linear divergences. Two-loop divergences cannot yet be handled, since 
we do not know what their precise form is. It may be interesting to study them, 
not only for the introduction of counterterms, but also to see if they have the 
same nice properties as the one-loop divergences (classical HTL's), such as gauge 
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invariance and a conserved energy for the effective theory. 
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A Hot, classical Feynman rules 
A.l Scalar fields 

In this appendix we discuss the construction of classical diagrams in perturbation 
theory, i.e. the classical Feynman rules at finite temperature, in scalar field theory. 
In order to do this, we start by repeating some necessary ingredients of the 
approach that was introduced in 0. For definiteness we use a massive scalar 
field with mass m, and interaction 

Sint = - J d^'x (^^g<P' + ^X<P'y (A.l) 

In classical perturbation theory, two type of two-point functions appear. A 
perturbative solution of the equations of motion, 

(9,^-V^ + m^)0 = ^, (A.2) 

is constructed with the free retarded Green function Gq{x), as 

<f>{x) = 0o(x) + / rfV G^{x - X )^^ + . . . , (A.3) 

where (poi^) is the solution of the unperturbed problem with some arbitrary initial 
conditions. Note that this yields expressions in which only 0o's are left over. The 
other two-point function specifies how to treat expectation values of 0o's- The 
free thermal propagator So{x) carries the thermal information and is defined by 

So{x - X') = (00 (a;) 00 (a;') )cl, free- (A.4) 

The brackets denote the averaging over the initial conditions weighted with the 
Boltzmann weight, for the unperturbed case In momentum space, the intro- 
duced two-point functions read 

^oiK) = ^2 _ + ,^)2 = E ^k^ + ,e + su^' ^^-^^ 
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So{K) = n,i{ky{k')27i6{kl -u;i)=Yl n,i{su^)-^27is6{k'> - su^), (A.6) 

s=± 

T 



nci{k') = j-^, cuk = Vk2 + m2, e(A;O) = ^(A;°)-0(-A;°). (A.7) 

The (free) retarded and (free) thermal two-point function are related by the 
classical KMS condition |^ 

tSoiK) = n,,{k^) [G^{K) - G^{K)\ , (A.8) 

where Gq{K) = Gq[—K) is the free advanced Green function. Finally, classi- 
cal loop integrals containing these two-point functions arise from the spacetime 
integral (s) in ( |A.3|) . 



Explicitly solving the equations of motion perturbatively and making all pos- 
sible contractions to find all possible diagrams becomes rather cumbersome at 
higher order in the coupling constants. Therefore we discuss in the remainder of 
this appendix a set of rules which are based on the underlying quantum pertur- 
bative approach. 

The imaginary-time or Matsubara formalism, does not lead to a close connec- 
tion with the classical approximation as described above at intermediate stages of 
a calculation. However, a useful observation, obtained using the imaginary-time 
formalism, is explained in section |3.1| . After performing the Matsubara sums, 
every diagram has a term which has as many distribution functions as loops. 
Hence, in the classical limit h ^ these terms remain. Other terms have less 
distribution functions and go to zero. This ensures us that the classical limit is 
non-trivial and exists. 

A formalism which lies closer to the classical perturbative approach is a vari- 
ation on the real-time formulation of finite temperature field theory, and uses 
the closed time path (CTP) method [Q. As is well-known, the CTP method 



involves a contour in the complex-time plane that consists of two branches, the 
upper branch and the lower branch C„ that runs back in time. This leads 
to a doubling of the fields, and they are denoted as 0+, 0_ to indicate on which 
branch they live. The propagator takes a matrix form, 

^, A ( G++ix-x') G+-ix-x')\ 

where the different superscripts specify the possible positions on and orderings 
along the contour. A convenient variation is based on the Keldysh formalism, 
and is the 'center-of-mass/relative' coordinates version. It uses a change of basis 
from 0+ „ to 01^2 

/ 01 \ / (0+ + 0_)/2 \ 

[h [ 0+ - 0- r 



(A.IO) 
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K 

(a) (b) 
Figure 5: Propagators, (a) G§{K) = G^{-K), (b) iFo{K). 





(a) 



(b) 



(c) 



(d) 



Figure 6: Vertices, (a) 5fi'0?02, (b) ^fi-^^, (c) h\(t)\<p2, (d) tMi4>1- 



such that the (free) matrix propagator takes the form 
G(o)(x - x') - 







(A.ll) 



Here the free retarded and advanced Green function are given in momentum space 
by the (classical) expression ( |A.5| ), and the quantum thermal two-point function 
in momentum space reads 



1, 1 



Z ZUJ- 



(A.12) 



which is of course the quantum version of (|A.6| ). Again the (free) retarded and 
thermal two-point functions are related by the KMS condition 



iFoiK) = n{k^) G^{K) - G^{K) 



(A.13) 



Feynman rules appear when also the interaction part along the closed time 



path contour is written in terms of the 0i^2 fields |^ 

5'int = 



d'x (^-gfy2 + \,g<t'l + + ^A^i^^) . (A.14) 

The rules are presented pictorially in figs. ^ and |. The 0i field is denoted with 
a full line and the <p2 field with a dashed line. For the retarded and advanced 
Green functions, it is necessary to specify the direction of the momentum fiow 
through the propagator, and this is indicated with the arrow. 
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(a) (b) 

Figure 7: (a) Retarded self-energy, (b) generalized retarded n-point vertex func- 
tions. 

The retarded self energy and the so-called generalized retarded n-point vertex 
functions |^ have one dashed 'leg' and n — 1 full 'legs'. These are shown in fig. 



1^. The arrows denote again the momentum flow of the external momenta. 

We now discuss the ^ — > limit of these real-time quantum Feynman rules. 
This limit affects the diagrams in two ways. The first one is obvious, the thermal 
propagator Fq has to be replaced by 5*0. The second change leads to a drastic 
simplification: only the vertices (a) and (c) in fig. ^ contribute in the classical 
limit, and the two other vertices (b) and (d) can be neglected. This can be seen 
as follows: vertices (b) and (d) can only appear in a diagram with retarded (or 
advanced) Green functions attached to the three dashed legs. After attaching 
these Green functions, the resulting outer lines (which either still have to be 
attached to another vertex or are external lines) are always full lines. However, 
such a configuration can be constructed as well with vertices (a) and (c): these 
vertices have two full legs where (b) resp. (d) have two dashed legs. By attaching 
two thermal two-point functions on these legs, the external lines are full as well, 
and the vertices can be part of a diagram in exactly the same manner. But a 
classical thermal two-point function is proportional to 1/h. Diagrams with vertex 
(a) or (c) have two more thermal two-point functions than the corresponding 
diagrams with vertex (b) or (d). Hence, the first class of diagrams is relatively 
stronger in the classical limit with respect to the second class by a factor 1/^^.0 
In other words, vertices (b) and (d) will be 0{h'^) suppressed with respect to 
vertices (a) and (c). 

We propose that classical Feynman rules follow from the quantum ones by 
taking h to zero, which results in the following (simple) rules: 

1. Draw all diagrams as in the quantum case, but use only vertices (a) and 
(c). 

2. Replace the thermal propagator Fq by its classical counterpart 5*0. 



^Negative powers of Ti will of course be canceled by positive powers coming from loop count- 
ing. 
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An explicit check of these rules (by a comparison with the results obtained by 
perturbatively solving the equations of motion and averaging over the initial 
conditions) can be found for the case of A(^^-theory for the two-point function up 
to two loops and the four-point function to one loop in 0. 

A. 2 Gauge invariant cut-ofF in the classical theory 

We argue that in classical gauge theories it is possible to introduce a (continuum) 
momentum cut-off without breaking gauge invariance. The basic ingredient is the 
result of Landshoff and Rebhan that in general linear gauges it is possible to 



formulate a (quantum) real-time theory in which only the two physical degrees 
of freedom of the gauge field acquire a thermal part. This implies that a change 
in the distribution function 

n{k'^) ^n{k^)f{k/A), (A.15) 

with / some function, does not break gauge invariance. Introducing a cut-off in 
this way will not affect the Slavnov- Taylor identities. This has been employed in a 
Wilson renormalization group approach to hot (quantum) SU(iV) gauge theories 

If we take the classical limit of ( |A.15| ) and choose / as the step function, we 



get 

n,ik')-^n,,{ky{A~k), (A.16) 



which as ( |A.15| ) does not break gauge invariance. It is for instance straightforward 



to check that the HTL's calculated with distribution function ( A.16|) satisfy the 



same abelian-like Ward identities as usual. Finally we should remark that the 
regularization ( |A.16|) is sufficient to render the theory ultraviolet finite, since each 
loop introduces one distribution function.^ 



B Classical one- loop SU(A^) self-energy: explicit 
calculation 

We present in this appendix the calculation of the classical self-energy in SU(A^) 
gauge theory, in particular the Iljj part, in the Feynman gauge. The starting point 
is given by (|2.5| ) in the main text. After changing variables from k — > — k — p in 
the part that is proportional to nd^tOp+k) , we find 

nti(^) = ^Vivnei(P), (B.i) 



In the quantum theory the cut-ofF in ( [A.15 ) acts only on thermal fluctuations. A zero- 



temperature regularization and renormalization is still necessary to avoid divergences coming 
from the zero-temperature quantum fluctuations. 
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with 

nei(P) 



6 + 



An 



(p°+^k)2 



p+k 



(B.2) 

and An = Ak"^ + 4k ■ p + — Gp^. We have combined h with nci(cc'k), which is 
an ^-independent combination. 

The angular integrations can be performed, and 



nci(/,p)= dk 



- / dk 



khnci{k) 



1 + ^ln^ 



P 



khuciik) r — 4pq r 



"I" 



k_ 
2^ 



L^{k)-L.{k) 



(B.3) 



L_(A;) +4p° 



Motivated by Weldon we used here the notation 



p± = -(p°±p), L±(A;)=ln 



A; ± 



(B.4) 



2'" "~k±p. 

The result ([B.3|) agrees with the expression obtained by Weldon in the appendix 
of P^l, except of course that the distribution function is classical in our case. 

The remaining radial integral in the first line of ( p.3|) is linearly divergent. For 
the first two terms this is obvious, and for the third term one can use L^{k) — 
L-{k) = 2p/k + 0{k~^). In fact, the divergence in this term cancels against 
the first term. The integrals in the second line are convergent. To regulate the 
divergences, we use the distribution function with a momentum cut-off hnc\{k) = 
T/k 6{A — k). The final result requires the evaluation of four integrals, which 
read (recall that p^ contains a small positive imaginary part) 



J dkkhriciik) = TA, 



dkk^hridik) L+{k) - L^{k) 



dkhriciik) L^{k) — L_{k) 
dk khriciik) \L+{k) + L_{k) 



T ( 2pA + -iripp 



Ttci In 



P- 
-Tmp. 



(B.5) 
(B.6) 
(B.7) 
(B.8) 



The second and fourth integral are straightforward using partial integration, and 
the third one can be performed by complex contour integration while being careful 
around k = 0. Note that these integrals are much simpler than in the quantum 
case, because of the simple k dependence of the classical distribution function. 

Putting all the results together, we find for the classical one-loop retarded self 
energy 

"^^ln^ + -^ ^-0 3p2-4p2 

tt'^ p p_ ^ 



nf 



which is presented in (P?^). 



ip 



2p 



-o-i P+ 
-I m — 

P- 



(B.9) 
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C Two loop naively linear divergent contribu- 
tions 



C.l Diagram b 



In this appendix we give the results for the naively linearly divergent contributions 
to the classical two-loop self-energy. We start with the classical limit of the self- 
energy diagram (b) in fig. |I|, presented in (^.51), and use the shorthand notation 
of (|3.6| ). There are three naively linearly divergent contributions and we shall 
denote these with (bl), (b2), and (b3). 

We start with contribution (bl), obtained by taking S3 = s',S2 = —s and 
setting the external p^, p to zero in the energy denominators with three loop- 
energies. We then find 



n 



(bl). 



Sk' 



-si 



(27r)3 ^ ''^ 2^uju'uJiUJ2UJ^ + s'{u;' + Us) 



-S'UJ' ~ sou + SiUJi 



S'uj' — SUJ + SiUJi 



■ (C.l) 



The difference between distribution functions [nc\{sujp_\^ — nc\{suj\.)] reduces the 
degree of divergence by one compared to the naive estimate, which is from linear 
to logarithmic. Note that the other difference between distribution functions 
[?2ci(sic<Jk-k') — '"■ci('S'^^p-k')]5 does not reduce the degree of divergence any further, 
since k is not a (small) external momentum, but is integrated over. 

A similar contribution is obtained by taking S2 = s and S3 = —s' and again 
setting p = in the same energy denominators. We obtain 



n(^!i(p) 



(t^k r d^k 



(2vr) 



1 



s'uj' + SUJ + SiUJi 



7^^ 2^ujuj'ujiUj2UJ2, p^ + s'{uj' — UJ^) 

(C.2) 



nc\{siuji) + nc\{suJ2) ^ rid^SiUJi) - rid^suj) 



s'uj' — SUJ + SiUJi 



Again a difference between distribution functions appears that reduces the degree 
of divergence to a logarithmic one. 

The third naively linearly divergent contribution to consider is of a different 
type. It is obtained from the classical limit of ( |3.5| ) by setting s = —S2 and 
s' = —S3 and taking the linear term in p^,p in an expansion of the energy 
denominator with uji = . The zeroth order term in this expansion gives rise 
to a naively quadratic divergence and was already discussed in the main text. 
The first-order term reads 



nSS(P) 



1 



d?k r <Pk 



{2Tif J (27r) 



2l'UJUj'uJ\UJ2UJ'i 



1 1 1 

+ S'{UJ' - iOs) p° + S{UJ - 002) {s'uj' - SUJ + SiUJif 

|s'(p • k') Ki(su;2) - ridisu)] [^^(s'cus) + n^iisiUi)] 

-s(p ■ k) [udis'uj^) - Udis'uj')] Ki(su;2) - ridisiUJi)] 
+pVi(sit^i)(Ki(st^) - nci{suj2)] + Ki(sV) - ncKs'cus)]) 

+pVi(sit^i) [nci{suj)nci{s'uJs) - nci{suj2)nci{s'uj')] |. (C.3) 

We emphasize again that the region of phase space where s'uj' — suj + siuji vanishes 
is excluded in this expansion. The first three terms between curly brackets all 
have a factor which is the difference between distribution functions. The fourth 
term is different, but also here the factor with distribution functions vanishes 
when the external momentum is taken to zero (i.e. when 002 uj,uj3 ^ u'). 
Hence this factor contributes a power instead of A~^, and it brings down the 
degree of divergence. We conclude that the degree of divergence is reduced from 
linear to logarithmic in contribution (b3) as well. 



C.2 Diagram c 

The final diagram that needs to be examined is diagram (c) in fig. 0. The quantum 
expression is 

^ ' ' J (2vr)3 J (27r)3 ^"'^^ ^^^^^ 2^u;u'iu,iU2 - su - 82^2 

'''n{suj) + l][n{s'uj') + l]n{siUJi) 



sou + s'uj' — SiUJi 

—n{suj)n{s'uj')[n{siUJi) + 1] 
+— ([n{s'uj') + l]n{siUJi)n{s2UJ2) - 

- S2UJ2 + s'uj' - SiUJi ^ 

n{s'uj')[n{siuji) + l][n{s2UJ2) + 1])}, (C.4) 

where in this case uji = (x'p-k-k' and we inserted {k)jj to indicate the two powers 
of momentum that come from the two three-point vertices. 

We take the classical limit of ( |(J.4D . The contribution with S2 = —s is naively 
linearly divergent, it reads 

r(c) 2+\2 f f d^^' (i\2 —s'Si 1 



(27r)3 J (27r)3 ' 2^ujuj'ujiuj2 - suj + su' 
— [n{suj) — n{suJ2)][n{siUJi) — n{s'uj')]. (C.5) 

SUJ -\- S UJ — SitUk+k' 

Again the first difference between distribution functions reduces the degree of 
divergence to a logarithmic one. 
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